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We consider a system of nonlinear spinor and a Bianchi type I gravitational fields in 
presence of viscous fluid. The nonlinear term in the spinor field Lagrangian is chosen to be 
XF , with X being a self-coupling constant and F being a function of the invariants / an J 
constructed from bilinear spinor forms S and P. Self-consistent solutions to the spinor and 
. BI gravitational field equations are obtained in terms of x, where x is the volume scale of 

BI universe. System of equations for x and e, where e is the energy of the viscous fluid, is 
deduced. This system is solved numerically for some special cases. 
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|>. i I. INTRODUCTION 

o 

The investigation of relativistic cosmological models usually has the energy momentum tensor 
of matter generated by a perfect fluid. To consider more realistic models, one must take into 
account the viscosity mechanisms, which have already attracted attention of many researchers. 
Misner suggested that strong dissipative due to the neutrino viscosity may considerably 

reduce the anisotropy of the black-body radiation. Viscosity mechanism in cosmology can explain 
the anomalously high entropy per baryon in the present universe [3, 4]. Bulk viscosity associated 
with the grand-unified-theory phase transition [[5D may lead to an inflationary scenario fiSH]. 

A uniform cosmological model filled with fluid which possesses pressure and second (bulk) 
viscosity was developed by Murphy [9]. The solutions that he found exhibit an interesting feature 
that the big bang type singularity appears in the infinite past. Exact solutions of the isotropic 



homogeneous cosmology for open, closed and flat universe have been found by Santos et al [10], 
with the bulk viscosity being a power function of energy density. 

The nature of cosmological solutions for homogeneous Bianchi type I (BI) model was inves- 
tigated by Belinsky and Khalatnikov 111 ill by taking into account a dissipative process due to vis- 
cosity. They showed that viscosity cannot remove the cosmological singularity but results in a 
qualitatively new behavior of the solutions near singularity. They found the remarkable property 
that during the time of the big bang matter is created by the gravitational field. BI solutions in case 
of stiff matter with a shear viscosity being the power function of of energy density were obtained 
by Banerjee [fl~2h , whereas BI models with bulk viscosity (rj) that is a power function of energy 
density e and when the universe is filled with stiff matter were studied by Huang II 1311 . The effect 
of bulk viscosity, with a time varying bulk viscous coefficient, on the evolution of isotropic FRW 
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models was investigated in the context of open thermodynamics system was studied by Desikan 
Iil4l1 . This study was further developed by Krori and Mukherjee [fl5ll for anisotropic Bianchi mod- 
els. Cosmological solutions with nonlinear bulk viscosity were obtained in II 1611 . Models with both 
shear and bulk viscosity were investigated in [17, Till . 

Though Murphy [9] claimed that the introduction of bulk viscosity can avoid the initial sin- 
gularity at finite past, results obtained in |fl9ll show that, it is, in general, not valid, since for 
some cases big bang singularity occurs in finite past. To eliminate the initial singularities a self- 
consistent system of nonlinear spinor and BI gravitational field was considered by us in a series 
of papers [|20l l2ll I22L 231. For some cases we were able to find field (both matter and gravita- 
tional) configurations those were always regular. In the papers mentioned above we considered 
the system of interacting nonlinear spinor and/or scalar fields in a BI universe filled with perfect 
fluid. We also study the above system in presence of cosmological constant A (both constant 
and time varying 12311 '). A nonlinear spinor field, suggested by the symmetric coupling between 
nucleons, muons, and leptons, has been investigated by Finkelstein et. al. ll24ll in the classical ap- 
proximation. Although the existence of spin- 1/2 fermion is both theoretically and experimentally 
undisputed, these are described by quantum spinor fields. Possible justifications for the existence 
of classical spinors has been addressed in Ii25ll . In view of what has been mentioned above, it 
would be interesting to study the influence of viscous fluid to the system of material (say spinor 
and/or scalar) and BI gravitational fields in presence of a A-term as well. In a recent paper we stud- 
ied the Bianchi type-I universe filled with viscous fluid in presence of a A term [26]. This study 
was further developed in [26] where we present qualitative analysis of the corresponding system 
of equations. Finally in [26] we introduced spinor field into the system and solved the system for 
some special choice of viscosity. The purpose of this paper is to further developed those results 
for more general cases and give some numerical results. It should be noted the in the process there 
occurs a very rich system of equations for volume scale, Hubble constant and energy density. The 
qualitative analysis of this system is under active study and we plan to present those results soon. 



II. DERIVATION OF BASIC EQUATIONS 

In this section we derive the fundamental equations for the interacting spinor, scalar and gravi- 
tational fields from the action and write their solutions in term of the volume scale T defined bellow 
(12.161) . We also derive the equation for % which plays the central role here. 

We consider a system of nonlinear spinor, scalar and BI gravitational field in presence of perfect 
fluid given by the action 

y(g; ¥ ,xj/)= I ^^dQ. (2.1) 

with 

=Sf — =$^g ~t~ =^sp =2m- (2.2) 

The gravitational part of the Lagrangian (12.21) is given by a Bianchi type I (BI hereafter) space- 
time, whereas Jz? sp describes the spinor field lagrangian and Jzf m stands for the lagrangian density 
of viscous fluid. 
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A. Material field Lagrangian 

For a spinor field i/a, symmetry between yr and \[r appears to demand that one should choose 
the symmetrized Lagrangian [29[]. Keep it in mind we choose the spinor field Lagrangian as 



sp 



my/yr + XF, 



(2.3) 



Here m is the spinor mass, A is the self-coupling constant and F — F(I,J) with I = S 2 = (yrys) 2 
and J = P 2 = [iyry^yr) 2 . According to the Pauli-Fierz theorem ||3~otl among the five invariants only 
/ and J are independent as all other can be expressed by them: Iy = = I + J and Iq = I — J. 
Therefore, the choice F = F(I,J), describes the nonlinearity in the most general of its form [21]. 
Note that setting A = in (12.31 ) we come to the case with linear spinor field. 



B. The gravitational field 

As a gravitational field we consider the Bianchi type I (BI) cosmological model. It is the sim- 
plest model of anisotropic universe that describes a homogeneous and spatially flat space-time 
and if filled with perfect fluid with the equation of state p = £e, £ < 1, it eventually evolves 
into a FRW universe rf3~fll . The isotropy of present-day universe makes BI model a prime candi- 
date for studying the possible effects of an anisotropy in the early universe on modern-day data 
observations. In view of what has been mentioned above we choose the gravitational part of the 
Lagrangian (12.21) in the form 

= (2-4) 

where R is the scalar curvature, K = SnG being the Einstein's gravitational constant. The gravita- 
tional field in our case is given by a Bianchi type I (BI) metric 

ds 2 = dt 2 -a 2 dx 2 -b 2 dy 2 -c 2 dz 2 1 (2.5) 
with a, b, c being the functions of time t only. Here the speed of light is taken to be unity. 



C. Field equations 

Let us now write the field equations corresponding to the action (12 .11 ). 
Variation of (12.11) with respect to spinor field yf[yf) gives spinor field equations 

*V i V M v/-mv/+^v/ + g ? ?'yV = 0, (2.6a) 
N^xj/y^ +m\j/- ^xjf-^ixj/y 5 = 0, (2.6b) 

where we denote 

dF dF 
^ = 2A5^-, <S = 2'kP%-. 
ol oJ 

Variation of (12.11) with respect to metric tensor g^ v gives the Einstein's equations which in 
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account of the A-term for the BI space-time (12.51) can be rewritten as 

b c b c i 

-+-+--= kT{ + A, (2.7a) 
b c be 

eded -> 

- + - + — = KTi + A, (2.7b) 
c a c a 

d b db ^ „„ „ N 

- + u + -u = kT 3+ A i (2.7c) 
a b ab 

db be cd n ,„ „ « 

-u + i:- + — = kT o+ A ^ ( 2 -7d) 
ao be c a 

where over dot means differentiation with respect to t and Ty is the energy-momentum tensor of 
the material field given by 

T^ = T S ^ + T m y (2.8) 

Here r sp v ^ is the energy-momentum tensor of the spinor field which with regard to (12.61) has the 
form 

= ^(^y^vW+wyv^^w-^iiWyvW-^vwyiiW^ (2.9) 

r^„ is the energy-momentum tensor of a viscous fluid having the form 

TZfi = ( e + P) u ii uV ~ P ^ + V8 V ^ \ u ii$ + u p-.fi ~ u^up-a - upu"u^ a ], (2.10) 

where 

P=p-tt-\r))u%. (2.11) 

Here £ is the energy density, p - pressure, r\ and £ are the coefficients of shear and bulk viscosity, 
respectively. In a comoving system of reference such that = (1, 0, 0, 0) we have 

r« = e, (2.12a) 

7^ = "/ + 2tA (2.12b) 

Til = -P' + 21^ (2.12c) 

7^3 = ~p'+2ri~. (2.12d) 



In the Eqs. (12.61) and (12.91) is the covariant derivatives acting on a spinor field as ||32l. 13311 

v MV = ^r-r M v, v M v/=^ + v/r M , (2.13) 

where are the Fock-Ivanenko spinor connection coefficients defined by 

r M = ^(rj^-fy/a). (2.14) 
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For the metric (12.51) one has the following components of the spinor connection coefficients 

To = 0, ^ = -a^yV, T 2 = l -b{t)ff, T 3 = l -c(t)ff. (2.15) 

The Dirac matrices y" (x) of curved space-time are connected with those of Minkowski one as 
follows: 

f = f, y x = t/a, f = f/b, r i = f/c 

with 

f-a°,\ MVo'V ^Mv; 



where a, are the Pauli matrices: 



l 1 2 -i\ 3 / 1 

° = V 1 J ' °-{io J' a = (,0-1 

Note that the y and the o matrices obey the following properties: 

ff + ff = 2r 1 i j, ij = 0,1,2,3 
?? + ?? = 0, (f) 2 = I, i = 0,l,2,3 
cWcr* = % + ie jk iO l , j, k, I = 1 , 2, 3 

where r/,/ = {1,-1,-1,-1} is the diagonal matrix, 5^ is the Kronekar symbol and e ;W is the 
totally antisymmetric matrix with £123 = +1. 

We study the space-independent solutions to the spinor field equations (12.61) so that 1// = y/(t). 
Here we define 

T = abc=y/ = g (2.16) 
The spinor field equation (|2.6al) in account of (12.131 ) and (12.151 ) takes the form 

if\ -=r + — J Y-m\j/ + @\j/ + &iy 5 \i/ = 0. (2.17) 
\ot 2t / 

Setting Vj(t) = y/xy/jit), j = 1,2,3,4, from (12.171) one deduces the following system of equa- 
tions: 

Vi + i{m-@)Vi-&V3 = 0, (2.18a) 

V 2 + i(m-3>)V2-yV 4 = 0, (2.18b) 

V 3 -i(m-gi)V 3 +&V 1 = 0, (2.18c) 

V4-i{m-@)V 4 + &V 2 = 0. (2.18d) 

From (|2.6al) we also write the equations for the invariants 5, P and A = i//y 5 y°i// 

S -2&A = 0, (2.19a) 
P -2(m-@)A = 0, (2.19b) 
A + 2(m-^)/ 3 o + 2^5 = 0, (2.19c) 

where So = zS, Pq = TP, and Aq = xA. The Eq. (12.191 ) leads to the following relation 

S 2 + P 2 +A 2 = C 2 /t 2 , C 2 = const. (2.20) 
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Giving the concrete form of F from (12.181) one writes the components of the spinor functions 
in explicitly and using the solutions obtained one can write the components of spinor current: 

f = ¥f¥- (2-21) 

The component j° 

f = - [v*vi + v^v 2 + v;v 3 + y; v 4 ] , (2.22) 

defines the charge density of spinor field that has the following chronometric-invariant form 

P = Uo-f) 1/2 . (2-23) 
The total charge of spinor field is defined as 

00 

Q = J P \f^gdxdydz = pxf , (2.24) 

— 00 

where V is the volume. From the spin tensor 

one finds chronometric invariant spin tensor 

4f =(%o^°) 1/2 , (2.26) 
and the projection of the spin vector on k axis 

& Q yf^gdxdydz = S%°xV. (2.27) 



Let us now solve the Einstein equations. To do it, we first write the expressions for the compo- 
nents of the energy-momentum tensor explicitly: 

(2.28a) 
(2.28b) 





= mS-XF + £ = 


J ) 




7? 


= @S + yP-XF 


-p' + 2r\- 
a 


= Tl+2r]-, 
a 


Ti 


= 3>S + &P-XF 








= &S + &P-XF 


-p'+2r}- 

c 





(2.28c) 
(2.28d) 

In account of (I2.28I ) subtracting (|2.7a|) from (|2.7b| ), one finds the following relation between a and 
b: 

-=D 1 exp^ i y J. (2.29) 

Analogically, one finds 

b ( fe- 2K ^ dt dt\ c ( fe- 2K ^ dt dt\ 
-=D 2 expfx 2 y J, -=D 3 cxp\X 3 j . (2.30) 
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Here D\, Z>2, D 3l X\, X2,X 3 are integration constants, obeying 

J DiD 2 D 3 = l, X 1 +X 2 + X 3 = 0. 



In view of (12.311) from (|2.29l) and (12.301) we write the metric functions explicitly [21] 



a{t) = ( J D 1 /D 3 ) 1 / 3 T 1 / 3 exp 



Xi-X 3 r e - 2K ^ dt 



b(t) 



c it) 



(DfD 3 )" 1/3 T 1/3 exp 



(D 1 D^) 1 / 3 T 1/3 exp 



3 

2X 1 +X 3 
3 

Xi + 2X 3 f e 



dt 

no 

e -2Kjridt 
-2k f r\dt 

-dt 



-dt 



x(t) 



(2.31) 

(2.32a) 
(2.32b) 
(2.32c) 



As one sees from (I2.32al) . (I2.32bl) and (I2.32cl) . for x = t n with n > 1 the exponent tends to unity at 
large t, and the anisotropic model becomes isotropic one. 

Further we will investigate the existence of singularity (singular point) of the gravitational 
case, which can be done by investigating the invariant characteristics of the space-time. In general 
relativity these invariants are composed from the curvature tensor and the metric one. In a 4D 
Riemann space-time there are 14 independent invariants. Instead of analyzing all 14 invariants, 
one can confine this study only in 3, namely the scalar curvature I\ = R, I2 = R^y/lV, and the 

Kretschmann scalar 73 = R a ^^ v R a ^^ v ■ At any regular space-time point, these three invariants 
7i, h : h should be finite. One can easily verify that 



h 



1 



h 



1 



h 



1 



Thus we see that at any space-time point, where T = the invariants /j, I2, 1 3 , as well as the scalar 
and spinor fields become infinity, hence the space-time becomes singular at this point. 

In what follows, we write the equation for x and study it in details. 

Summation of Einstein equations (I2.7al) . (I2.7bl) . (I2.7cl) and (I2.7dl) multiplied by 3 gives 



x =?-k(t§ + T?)x + 3kt]x + 3Ax, 



which can be rearranged as 

3 



x - - K ^x = -K(mS + @S + <3P - 2XF + £ - p^jx + 3Ax. 



(2.33) 



(2.34) 



For the right-hand-side of (|2.34l) to be a function of x only, the solution to this equation is well- 
known [|34f| . 

On the other hand from Bianchi identity G^. v = one finds 



which in our case has the form 



(2.35) 



(2.36) 
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ho - (a 2 b 2 c 2 



This equation can be rewritten as 

f„»=-^-^ j+ 2n^ + F + ? j. (2.37) 
Recall that (12.191 ) gives 

(m-&)So-&P = 0. 
In view of that after a little manipulation from (12.371 ) we obtain 

£ + 1(0-^ + U)^ +477(k7 ° + A) = 0, (2.38) 

where 

co = e+ j p, (2.39) 

is the thermal function. For further purpose we would like to note that in absence of shear viscosity 
from Eqs. (12.331) and (12.371) one finds 

fcf ° = 3H 2 - A + Coo, Coo = const. (2.40) 

where in analogy with Hubble constant introduce the quantity H, such that 

* = * + ± + 6 - = 3H. (2.4i) 

x a b c 

Then (12.341) and (12.381) in account of (12.281) can be rewritten as 

H = l(3£H-a>)- (3H 2 - ?ce - -A) + 1 (mS + &S + &P- 2XF) , (2.42a) 

e = 3H(3^H-co)+4r](3H 2 ~K£-A)-4r]K(mS-?iF). (2 Alb) 

Thus, the metric functions are found explicitly in terms of x and viscosity. To write x and com- 
ponents of spinor field as well and scalar one we have to specify F in Jzf sp . In the next section we 
explicitly solve Eqs. (12.181) and (|2.42l) for some concrete value of F . 

The Eqs. (|2.42l) can be written in terms of dynamical scalar as well. For this purpose let us 
introduce the dynamical scalars such as the expansion and the shear scalar as usual 

= u^ o 2 = l -o llv o^ L \ (2.43) 

where 

<V = 2 ( m m;« p v +u V;a P") ~ ^OPiiv (2.44) 
Here P is the projection operator obeying 

P 2 = P. (2.45) 
For the space-time with signature (+, — , — , — ) it has the form 

Pnv =gixv-"M"v, Py = 8$ -u^u v . (2.46) 
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For the BI metric the dynamical scalar has the form 

d b 

- + T + 
abet 



a b c x 

6 = - + - + - = - (2.47) 



and 



a 



b 2 c 2 3 



In account of (12.321) one can also rewrite share scalar as 



2o 2 = 6( * ? + *f +X 'K- 4K ^ dl . (2.49) 



From (|2.7dl) now yields 

l -0 2 -a 2 = K 



mS-XF + e 



+ A (2.50) 



The Eqs. (|2.42l) now can be written in terms of and a as follows 

= ^.^e-0))-^-(mS-^S-^P)-3o 2 1 (2.51a) 

e = 0(^0 -co) +4r)<T 2 . (2.51b) 
Note that the Eqs. (12.511) without spinor and scalar field contributions coincide with the ones given 



in H2D 



III. SOME SPECIAL SOLUTIONS 

In this section we first solve the spinor field equations for some special choice of F, which will 
be given in terms of %. Thereafter, we will study the system (12.421) in details and give explicit 
solution for some special cases. 

A. Solutions to the spinor field equations 

As one sees, introduction of viscous fluid has no direct effect on the system of spinor field 
equations (12.181) . Viscous fluid has an implicit influence on the system through T. A detailed 
analysis of the system in question can be found in [|2l|l . Here we just write the final results. 

1. Case with F =F{I) 

Here we consider the case when the nonlinear spinor field is given by F = F(I) . As in the case 
with minimal coupling from (I2.19al) one finds 

S=—, Co = const. (3.1) 

T 



For components of spinor field we find 112 111 



v T 



(3.2) 
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with C, being the integration constants and are related to Co as Co = C 2 + C| — C 2 — C|. Here 
P=f(m-@)dt. 

For the components of the spin current from (12.211) we find 

f = I[ Cl +C 2 2 + C 3 2 + C 4 2 ], j 1 = ^[C 1 C 4 + C 2 C 3 ]cos(2/3), 
j 2 = A [C1C4 -C2C3] sin(2j3), / = ^ [C1C3 -C2C4] cos(2j8), 
whereas, for the projection of spin vectors on the X, Y and Z axis we find 

£23,0 _ C1C2+C3C4 31)0 _ q ^12,0 _ ^1 ~ ^2 + ^3 ~ ^4 

The total charge of the system in a volume Y in this case is 

Q = [C 2 + C 2 + C 2 +C 2 }y. (3.3) 

Thus, for T 7^ the components of spin current and the projection of spin vectors are singularity- 
free and the total charge of the system in a finite volume is always finite. Note that, setting A = 0, 
i.e., /3 = mt in the foregoing expressions one get the results for the linear spinor field. 

2. Case with F = F [J) 

Here we consider the case withF =F(J). In this case we assume the spinor field to be massless. 
Note that, in the unified nonlinear spinor theory of Heisenberg, the massive term remains absent, 
and according to Heisenberg, the particle mass should be obtained as a result of quantization of 



spinor prematter 113711 . In the nonlinear generalization of classical field equations, the massive 
term does not possess the significance that it possesses in the linear one, as it by no means defines 
total energy (or mass) of the nonlinear field system. Thus without losing the generality we can 
consider massless spinor field putting m = 0. Then from (|2.19b| ) one gets 

P = D /x, D = const. (3.4) 

In this case the spinor field components take the form 

Yi = 4= (D\e ia + iD 3 e- ia ) , yf 2 = -J= i D ^ G + iD^ ) , 

V X 



(3.5) 



y/3 = 4= (iDie ia +D 3 e- ia ) , i// 4 = 4= (* D 2^ +D A e- iG ) . 



The integration constants D, are connected to Do by Do = 2 (D 2 + D 2 ,— D 2 — D^). Here we set 
o = J^dt. 

For the components of the spin current from (12.211) we find 

f = l[ D \ +D l + Dl +D ^ : / = — [D 2 D 3 +D 1 D 4 ]cos(2a), 
x ax 

j 2 = — [D 2 D 3 -D1D4] sin(2cr), / = — [D1D3 -D 2 D 4 ] cos(2<x), 
bx cx 
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whereas, for the projection of spin vectors on the X, Y and Z axis we find 

s23 ,o_ 2(D l D 2 +D 3 D4) ^ S 3i, = s i2.o_ Dj-D 2 2+D 2 3-Dl 
hex 2abx 

We see that for any nontrivial x as in previous case the components of spin current and the pro- 
jection of spin vectors are singularity-free and the total charge of the system in a finite volume is 
always finite. 



B. Determination of x 



In this subsection we simultaneously solve the system of equations for x and £. Since setting 
m = in the equations for F = F(I) one comes to the case when F = F(J), we consider the case 
with F being the function of / only. Let F be the power function of S, i.e., F = S n . As it was 
established earlier, in this case S = Cq/x, or setting Co = 1 simply S = l/x. Evaluating S> in terms 
of x we then come to the following system of equations 



3k „ 3k (m X(n — 2) 
x = T ^ + T (- + -L_J + e _ p)T + 3 AT, 



4 A 1 



+ + -7])-3-4t] 



( m h . 



or in terms of H 



3Hx, 



1 , „ \ / o \ K (m X(n — 2)\ 

H = -(3^H-co)-(3H 2 -Ke-A) + -(- + ^ r ^). 



£ = 3H(3%H-Q))+4ri(3H 2 -K£-A)-4r}K 



m A 
X x n . 



(3.6a) 
(3.6b) 

(3.7a) 
(3.7b) 

(3.7c) 



Here r\ and % are the bulk and shear viscosity, respectively and they are both positively definite, 
i.e., 

T7 > 0, t, > 0. (3.8) 
They may be either constant or function of time or energy. We consider the case when 



ri=A£ a , %=B£P, 
with A and B being some positive quantities. For p we set as in perfect fluid, 

p = Ce, Ce(0,l]. 



(3.9) 



(3.10) 



Note that in this case £ ^ 0, since for dust pressure, hence temperature is zero, that results in 
vanishing viscosity. 

The system (13.71) without spinor field have been extensively studied in literature either partially 
ifiEtLl] or as a whole [11]. Here we try to solve the system (13.61) for some particular choice of 
parameters. 
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1. Case with bulk viscosity 

Let us first consider the case with bulk viscosity alone setting coefficient of shear viscosity 
77 = 0. We also demand the coefficient of bulk viscosity be inverse proportional to expansion, i.e., 

^q = ^H = C 2i C 2 = const. (3.11) 

Inserting r\ = 0, (|3.11|) and (13.101) into (|3.7cl) one finds 



1 



i + C 



C 3 = const. (3.12) 



Then from (I3.6al) we get the following equation for determining x: 

3k rC 2 I 3k(1-C)C 2 t 1+ £-C3 3kA(/z-2) . 

where g is the set of problem parameters. As one sees, the right hand side of the Eq. (13.131) is 
a function of x, hence can be solved in quadrature Q34{] . We solve the Eq. (|3.13l) numerically. It 
can be noted that the Eq. (13.131) can be viewed as one describing the motion of a single particle. 
Sometimes it is useful to plot the potential of the corresponding equation which in this case is 

W(q,x) = -lj&(q,x)dx. (3.14) 

The problem parameters are chosen as follows: K = 1, m = 1, A = 0.5, £ = 1/3, n = 4, C2 = 2 
and C3 = 1. Here we consider the cases with different A, namely with A = —2,0, 1, respectively. 
The initial value of x is taken to be a small one, whereas, the first derivative of x, i.e., t at that 
point of time is calculated from (12.401) . In Fig. [T]we have illustrated the potential corresponding to 
Eq. (13.131) . As one sees, independent to the sign of A we have the expanding mode of evolution, 
though a positive A accelerates the process, while the negative one decelerates. Corresponding 
behavior of x is given in Fig. [2] 



2. Case with bulk and shear viscosities 



Let us consider a more general case. Following |26(] we choose the shear viscosity being 
proportional to the expansion, namely, 

In absence of spinor field this assumption leads to 

3H 2 = K£+C 4 , C A = const. (3.16) 

It can be shown that the relation (13.161) in our case can be achieved only for massless spinor field 
with the nonlinear term being 

f = FqS 2{ - k - 1 ^ k . 

Equation for x in this case has the form 

Tf-0.5(l-C)t 2 -1.5^Tt-3[A-0.5(l-C)C 4 -AF T 2 ( 1 -^/' f ]T 2 = 0. (3.17) 
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FIG. 1: View of the potential corresponding to FIG. 2: Evolution of z depending on the signs of 

the different sign of the A term. the A term. 



In case of t, = const, and A = there exists several special solutions available in handbooks on 
differential equations. But for nonzero X we can investigate this equation only numerically. We 
consider the case when the bulk viscosity is given by a constant. Taking this into account for 
problem parameters we set £ = 1/3, ^ = 1, Fq = 1, X = 0.5 and C4 = 1. We study the role of A 
term. In doing this we consider the cases with positive, negative and trivial A. Since the nonlinear 
term in this case depends of fc, we also consider the cases with different K, namely with K > 1 and 
K < 1. In Figs. [3] and H] the evolution of % is illustrated for fc < 1 and fc > 1, respectively. In case 
of K < 1 we have non-periodic mode of evolution for all A, while for K > 1 a negative A gives a 
non-periodic mode of expansion. A non-negative A in this case gives an ever expanding mode of 
evolution. 



IV. CONCLUSION 

We consider a self consistent system of nonlinear spinor and gravitational fields within the 
framework of Bianchi type-I cosmological model filled with viscous fluid. The spinor filed non- 
linearity is taken to be some power law of the invariants of bilinear spinor forms. Solutions to the 
corresponding equations are given in terms of the volume scale of the BI space-time, i.e., in terms 
of T = abc. The system of equations for determining T, energy-density of the viscous fluid e and 
Hubble parameter H has been worked out. Exact solution to the aforementioned system has been 
given only for the case of bulk viscosity. As one sees from (12.421) or (12.511) . the system in question 
is a multi-parametric one and may have several solutions depending on the choice of the problem 
parameters. As one sees, solutions can be non-periodic independent to the sign of A term. Given 
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FIG. 3: Evolution of the universe with nontrivial FIG. 4: Evolution of the universe for different 

A term and K < 1. values of A term with K > 1. 



the richness of (13.61) we plan to give qualitative analysis of this system in near future. 
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